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1. Introduction 

The appearance of interfaces in statistical systems under particular conditions has always 
raised much interest in various fields of research, ranging from condensed matter to high 
energy physics. Recently, there have been remarkable theoretical and computational im- 
provements in the study of this phenomenon. 

Numerical investigations of interfaces in statistical systems (and of their analogue 
in Lattice Gauge Theories) have attained a great level of accuracy and reliability. In 
particular, the interface free energy in the Ising 3d model has been thoroughly studied by 
means of Monte Carlo simulations. Indeed, spin models provide a simple realization of 
interfaces since in their broken symmetry phase an interface separating coexisting vacua 
of different magnetization can be easily induced in the system by suitably choosing the 
boundary conditions. 

On the theoretical side, different effective models can be used to describe the behaviour 
of interfaces in 3d systems and in particular to evaluate their free energy. The most popular 
is the capillary wave model (CWM)[1, 2] which is based on the assumption of an action pro- 
portional to the area of the surface swept by the interface (for a review, sec for instance [3]). 
This model is tantamount to consider the interfaces as bosonic strings embedded in three 
dimensions, with a Nambu-Goto action [4, 5]. Although this effective string approach ne- 
glects the conformal anomaly which appears for target space dimensions d ^ 26, the effects 
of this approximation appear to be subleading [6] for large worldsheets. This description 
has exactly the same nature of the effective string description of certain observables, such 
as Polyakov and Wilson loops, in LGT (see for instance [7, 8] and references therein). 

The standard procedure to treat the effective Nambu-Goto string [9] in the last twenty 
years was to fix the so called "physical gauge" and expand perturbatively the partition 
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function around the classical solution (the surface of minimal area) in terms of the inverse 
product of the string tension a with the area of the surface. 

The increasing precision of recent numerical simulations allowed in the last few years 
to test this expansion beyond the leading order [10, 11]. Resorting to the first order formu- 
lation a la Polyakov of the NG model, it has been possible to re-sum the loop expansion 
obtained in the physical gauge- fixing, in particular in the case of Polyakov correlators [12] 
and of interfaces [13]. This treatment of the NG model yields a very good agreement with 
Monte Carlo results, at least for sufficiently large world-sheets; from this point of view, the 
comparison with MG data sets establishes a lower bound on the world-sheet area, below 
which one cannot neglect the effects of the conformal anomaly. 

Our paper lays in this context. Wc focus on the universal properties of periodic in- 
terfaces arising in general 2d models. Following the philosophy of the CWM we make the 
simple assumption that the weight of the interface is proportional to its length (which in- 
deed corresponds to a 2d version of the CWM) and treat at the quantum level this action 
by means of its first-order formulation. In this way we show in section 2 that the dominant 
term in the partition function acquires an universal form, proportional to to LmKi{mR), 
where m is the inverse of the correlation length, R is the lattice size in the direction of 
the interface, L the lattice size orthogonal to the interface and Ki is the modified Bessel 
function of the first order. 

In section 3 we consider an explicit statistical model, the 2d Ising model. Upon a 
suitable choice of boundary conditions, this system allows the formation of interfaces. The 
2d Ising model is under full analytic control, and we can derive directly the form of the 
dominant term in the interface partition function. We find agreement with the universal 
expression predicted by our 2d CWM. 

In section 4, we consider the exact expression of the free energy for interfaces in 3d 
obtained in [12], and we perform a dimensional reduction in one of the directions along 
the interface. We retrieve in this way the 2d expression proportional to Ki{mR) which we 
propose in this paper. This dimensional reduction is the exact analogue for the interface 
boundary conditions of the dimensional reduction from Polyakov loop correlators to spin- 
spin correlators studied in [14, 15]. A non trivial consistency test of this correspondence 
is that the relation which links the 2d mass m with the 3d string tension a is the same in 
the case of interfaces and of Polyakov loops. 

Recently, a new set of high precision Monte Carlo data for the free energy of interfaces 
in the 3d Ising model became available [16]. These data include asymmetric geometries in 
which one of the sides of the interface becomes much smaller than the other; in this situation 
one expects that the dimensionally reduced expression can describe the data accurately. 
In the last section, we test this expectation comparing these data with the predictions of 
the full 3d Nambu-Goto treatment and of the 2d simple model described here. 

2. A simple model for 2d interfaces 

Consider a physical system, defined on a two-dimensional space with periodic boundary 
conditions in both directions, in which two different phases, separated by a one-dimensional 
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Figure 1: An interface (the heavy magenta line F) separating two phases (blue and red) on a 
two-dimensional torus. 

interface, may coexist (see Fig. 1). 

We will describe this situation in a very simple way, analogous to the so-called capillary 
wave model [1, 2] for two-dimensional interfaces in 3d systems. We make the assumption 
that the effective action is just given by the length of the interface T: 



where r parametrizes the curve F, the dots denotes r derivatives, (x)^ = XjX* and m is 
a constant^ with the dimension of an inverse length. This is just the analogue of the 
conjecture which, in three dimensions, leads to the CWM [1, 2]. 

The corresponding partition function is given by the functional integral 



and it should correspond to the interface free energy, in analogy to what happens in the 
stringy description of 3d interfaces [13] . 

The action eq. (2.1), as well known, admits a first-order reformulation: 




(2.1) 






(2.3) 



where e(r) is the einbein^. 



so that the partition function can be written as 




(2.4) 



The action eq. (2.3) is invariant under reparametrizations of T: 



6t = e 



5e = ee + ee . 



(2.5) 



^We can view m as the "mass" of a relativistic particle in an Euclidean target space (the torus) whose 
world-line is the interface. 

^Namely one has ds^ — e^(r)dr^ 
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We can gauge-fix this invariance by choosing e(r) = A = const, so that the action becomes 

S[x,e]^'^ + ^j\r{x)^. (2.6) 

The constant A represents the length of the path and is a Teichmiiller parameter: it cannot 
be changed using the reparametrizations in eq. (2.5). As such, it must be integrated over 
in the partition function, which now reads 

mX 



Here 



dXe'—Z^Zgh- (2.7) 

Z^ = y Dor exp y' dr(x)2| (2.8) 

and we took into account the Faddeev-Popov ghosts b, c for the chosen gauge-fixing of the 
invariances in eq. (2.5). For them we have, by standard treatment, 

Zgh = DbDc exp j-A^ drbdrc^ . (2.9) 
With the chosen periodic boundary conditions on r, it turns out that 

Zgh = ^ . (2.10) 

The partition function is simply that of a non-relativistic particle of mass /x = m/A, 
whose Hamiltonian is 

i^=g. (2.11) 

Since the particle moves on a two-dimensional torus of sides Li,L2, the components of the 
momentum are quantized: 

P' = ^ . (2.12) 

for i = 1,2, with G Z. The partition function of this system is 

In the second step above we performed a Poisson re-summation, which reorganizes the 
sum in contributions of sectors describing quantum fluctuations around classical solutions 
of wrapping numbers mj. 

To describe an interface such as the one depicted in Fig. 1, we set 

mi = , m2 = 1 , (2.14) 



^Notice that in the first-order formulation we have to take into account the constraints which follow 
from the variation of 5 w.r.t. the einbein e, which (in the chosen gauge) are simply (x)^ = A. 
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getting therefore (if we remember that /x = m/A) 

Z. = L,L,^exp{-^Li} . (2.15) 

Inserting eq. (2.10) and eq. (2.15) into the expression eq. (2.7) of the total partition 
function we obtain finally^ 

where Ki is the modified Bcsscl function of first order. 

In the next section, we will consider a concrete physical system where the interface 
free energy can be computed directly, namely the 2d Ising model, and we will find that our 
proposed general expression eq. (2.17) does indeed capture its dominant behaviour. 

3. Interfaces in the 2d Ising model 

Interfaces can be realized in the context of spin models, since under peculiar boundary 
conditions their broken symmetry phase allows separated but coexisting vacua of different 
magnetization. Being interested in interfaces in 2d systems, it is natural to consider the 
2d Ising model, where an explicit analytical treatment is viable. 

To evaluate the interface free energy of the 2d Ising model in a simple way, we can 
take advantage of the fact that the continuum limit of the Ising model is described by 
the Quantum Field Theory of a free fermionic field of mass m. This framework allows the 
calculation of the partition function with any type of boundary conditions. Being interested 
in estimating the interface free energy, we will impose periodic boundary conditions in one 
direction and antiperiodic ones in the other direction. Explicit expressions of various 
partition functions can be found in [17] and (using the Q function regularization) in [18]^. 
It's worth pointing out that the solution proposed in [17] is nothing else than the continuum 
limit of Kaufmann solution [19]. Such continuum limit was performed by Ferdinand and 
Fischer in [20] at the critical point and the result of [17] is essentially the extension of [20] 
outside the critical point. 

Using the notations of [17], eq.s (91)-(105), the partition function of a fermion on a 
rectangular torus of sizes Li and L2, which for notational simplicity we shall henceforth 
denote as L and R, is given by the so-callcd "massive fermion determinant" : 

L>a,^(m]L,i?) = e-'5'"^^'^'5('')/6-^ J| (1 - (5«e-^^"('")/^) . (3.1) 

neZ+/3 

Here a, /3 can take the values 0, 1/2 and label the boundary conditions in the L and R 
directions respectively: the value corresponds to periodic b.c.s, 1/2 to antiperiodic ones. 

''We use the integral 

^~«e— f ^.(|)-V,_.,2^i,). (2.16, 

^Notice however that in [18] there is a sign mistake which was later corrected in [17]. 
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^nW_./^2 + (?^)2. (3.2) 



Moreover, 5a = e and 

e 

~ R V" ' ' R 

The adimensional variable r = mR sets the scale of the theory. 

Depending on the boundary conditions, the coefficient cp{r) can assume the following 
values: 



-\ir) = -,Y.^-^K,{kr) (3.3) 



2 71-2 Z__^ ^ 

A;=l 



or 

oo 



fc=l 

The partition function of the Ising model with periodic boundary conditions in both 
directions can be written as follows: 

Zising(m|L,i?) = ^ (Di 1 +L>o_i +Di o-L>o,o) . (3.5) 
If the L direction is antiperiodic, the partition function becomes 

^ising,ap(m|L, R) = \ {Di_i_ + D^ i_ - Di Q - Do,o) . (3.6) 

This is the situation which generates an odd number of interfaces along the R direction. 
The interface free energy will be given by the following expression: 

e-^inte.ace ^ Z,^,,,,^, = . (3.7) 

Ising 

Notice, as a side remark, that taking the limit^ r ^ 0, one moves to the critical point 
and all the expressions written above flow toward the conformal invariant ones and agree 
with the standard CFT results and with those reported in the pioneering work of Ferdinand 
and Fisher [20]. 

We are now interested in the large mL and mR limit of cq.(3.7). The large mL limit 
allows us to neglect the infinite product in eq.(3.1) while in the large mR we approximate 

fir 

ci(r)~co(r)~ ^ifi(r) . (3.8) 
The partition functions in eq.(3.5) and (3.6) simplify and we end up with 

-Rising ^ 1 . 



^Ising,ap ^ 



ttL / , , , N (3.9) 
ci{r) + co{r)j , 



from which we find 



ttT 9 r r 9 r m 

^interface ~ 77^(^1 (r) + co(r)) ~ — ^i^i(r) = -^K,{mR) . (3.10) 

By comparing this result to eq. (2.17), we see that the behaviour of the interface free 
energy in the 2d Ising model for large scales is perfectly captured by our simple effective 
model. 

®In this regime, one must keep in mind that the parameter t of [20] corresponds to r/2. 
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4. Dimensional reduction of the Nambu-Goto effective description of 3d 
interfaces 

In this section we shall discuss the behaviour of the Nambu-Goto effective string description 
of fluctuating interfaces when dimensional reduction occurs along one of the two lattice sizes 
which define the interface. 

As we discussed above, interfaces can be realized in spin models. In this realization, 
the physical ideas underlying the dimensional reduction process are perhaps best discussed. 
Focusing on the 3d Ising model, let us recall that it is mapped by duality into the 3d Ising 
gauge model. Under this mapping, the interface free energy becomes a close relative of the 
Wilson loop expectation value of the 3d gauge Ising model, the only difference being in the 
boundary conditions of the two observable: fixed in the Wilson loop case and periodic in 
the interface case. 

In lattice gauge theories, dimensional reduction has a very important physical meaning. 
The size of the lattice direction which is reduced is interpreted as the inverse temperature of 
the system and the point in which dimensional reduction occurs coincides with the decon- 
finement phase transition. The behaviour of the model in the vicinity of the deconfinement 
transition is well described, according to the Svetitsky-Yaffe conjecture [21], by the 2d Ising 
spin model. More generally, for other gauge groups with continuous deconfinement phase 
transitions, it is described by the 2d spin model with global symmetry group the center of 
the original gauge group. 

As a consequence of this observation, we expect that the interface free energy of the 3d 
Ising model should smoothly map into the 2d one when the dimensional reduction scale is 
approached from above, while it should vanish below it, since in the dual model this is the 
deconfined phase. Any reliable effective description of the free energy should therefore be 
compatible with these two requirements. This is a rather non trivial test for the effective 
string models which are expected to effectively describe both the Wilson and Polyakov loop 
expectation value and the interface free energy. 

In this section, we start from the expression for the interface free energy in 3d obtained 
in [13] and consider its behaviour under dimensional reduction; we find that it does indeed 
smoothly reduce to the 2d expression given here in eq. (2.17). This represents a test of 
our effective approach to the interface free energy both in three and in two dimensions. 

In the next section, we will refine this test by comparing both the 3d and the 2d the- 
oretical expressions to Monte Carlo data for asymmetric interfaces: the more asymmetric 
are the interfaces, the more we expect the prediction of our 3d and 2d expression to become 
compatible and reliable. 

Let us consider a three-dimensional torus of sides i^i,2,3, and assume that the interface 
lies in the plane orthogonal to Li. Our starting point is the expression for the 3d interface 
partition function in the Nambu-Goto approximation evaluated in [13] (eq.s (2.23, 2.24)): 




oo 



(4.1) 
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where 

A = L3L2 , u=-^ (4.2) 
^3 

and 



is a function of the occupation numbers and k' . The coefficients and c^/ are the number 
of partitions of k and k' . 

The dimensional reduction of eq. (4.1) can be most easily discussed by writing the 
interface free energy as a function of L3. Following the notations of [13], let us introduce 
the "energy levels" Ek^k' = crLsS. With this choice the argument of the Bessel functions 
in eq.(4.1) becomes L2-Efc,jt' and we have: 



Ek,k' = y cx^Li + Ana{k + k' - ^) + ■ (4.4) 

From this expression it is evident that, as the ratio increases, the separation between 

consecutive levels in the spectrum becomes larger and larger and, in particular, the gap 
between the lowest state k = k' = and the second one increases. 

In the dimensional reduction limit L3 ^ L2 we can thus truncate the sum in eq.(4.1) 
to the first term. In this limit, the argument of the Bessel function in eq.(4.1) becomes 
large, allowing us to use the asymptotic expansion 



1 + 0{1]}, (4.5) 



from which we see that higher states are exponentially suppressed with respect to the 
k = k' = one. Setting A; = A;' = 0, we get 



where now 



"«=,l-^- (1.7) 



3a A" V 3a Ll 



Equation (4.6) can therefore be rewritten as: 



^(0,0) 



1 



If we now define 



we find 



meff = aLa J 1 - (4.9) 



1 

^(0,0) = [l) ' ^1 "^effi^l (meffL2) , (4.10) 
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which is in perfect agreement with the partition function given by (2.17) a part from the 
normahzation constant. 

It is clear from eq. (4.7) that the above discussion is consistent only for 



In the string language, this bound represents the tachyonic singularity of the bosonic string: 
for L3 < the lowest state has a negative energy. In the 3d Ising model framework, 

the critical value 



at which the mass of the lowest state vanishes can be considered as the effective string 
prediction for the dimensional reduction scale which we were looking for the dimensional 
reduction scale which we were looking for. We thus see that this scale is naturally embedded 
in the NG formulation of the interface free energy. The same argument in the dual case 
of the Polyakov loop correlator leads to the identification of this scale with deconfinement 
temperature [22]. The value obtained in this way for the dimensional reduction scale 
(or, equivalently, for the deconfinement temperature in the dual model) is in rather good 
agreement with the Monte Carlo estimates. For the 3d gauge Ising model the disagreement 
is only of about 20%, as the Monte Carlo estimate turns out to be y/aLc = 0.8186(16) [23] 
instead of 1.023...) and it further decreases if one looks at SU {N) pure Yang Mills theories. 

It is important to stress that the mapping between 3d and 2d obscrvables which emerges 
from this comparison is exactly the same which is found considering the finite temperature 
behaviour of the Polyakov loop correlators in the vicinity of the deconfinement transition 
in the 3d gauge Ising model. In particular, the relation between the mass scale of the 
2d model and the string tension (in the present case the interface tension) of the 3d one 
reported in eq. (4.9) is the same which one finds in the Polyakov loop case. This is a rather 
non trivial consistency check of the whole dimensional reduction picture. 

As a last remark let us stress that the smooth flow of the 3d Nambu-Goto result 
towards the 2d ones under dimensional reduction has various important implications on 
our understanding of the Nambu-Goto effective action. Recall that the treatment leading 
to eq. (4.1) neglected the conformal anomaly, i.e., equivalently, did not take into account 
the Liouville field which does not decouple in the Polyakov formulation of the string when 
d / 26. Most probably the smooth fiow toward the 2d result occurs exactly because we 
neglected the Liouville field contribution. It is thus interesting to compare simultaneously 
the Nambu-Goto expression (without Liouville field) , the 2d result presented here and the 
output of the Montecarlo simulations for the free energy of interfaces on lattices with a 
value of L3 small enough to make the 2d approximation a reliable one. We shall devote 
the next section to this comparison. 

5. Comparison with the numerical data 

Following the above discussion we compare in this section the expressions for the 2d inter- 
face free energy given in eq. (2.17) and for the 3d free energy [13], reported in eq. (4.1), 




(4.11) 




(4.12) 
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64 
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10.6278 
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11.3138(25) 


11.3209 


11.3214 


11.3209 


11.3548 


11.3213 



Table 1: Numerical results of the free energy in the 3d Ising model at /3 = 0.223101 are compared 
to various theoretical estimates. See the main text for detailed explanations. 



^3 


L2 


U 


num 


F 


-^(0,0) 


-Fist 


-^1— loop 


-F2-I00P 


24 


64 


96 


18.4131(26) 


18.4121 


18.4121 
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96 


21.2414(27) 


21.2450 
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24.0306 


24.0306 


24.0306 


24.0497 


24.0312 


36 


64 


96 


26.7859(28) 


26.7873 


26.7873 


26.7873 


26.8017 


26.7875 


40 


64 


96 


29.5271(30) 


29.5250 


29.5251 


29.5250 


29.5365 


29.5251 


44 


64 


96 


32.2449(31) 


32.2498 


32.2500 


32.2498 


32.2594 


32.2498 


48 


64 


96 


34.9623(33) 


34.9653 


34.9657 


34.9653 


34.9736 


34.9653 



Table 2: Numerical results of the free energy in the 3d Ising model at /? = 0.226102 are compared 
to various theoretical estimates. See the main text for detailed explanations. 

with a set of high-precision Monte Carlo data [16] obtained in a 3d Ising model in which 
interfaces are created by imposing antiperiodic boundary conditions in one of the three 

lattice directions. 

In particular we isolated two sets of data. The first one is reported in Tab. 1 and 
displays the values of the free energies obtained in the 3d Ising model at /3 = 0.223101. 
The interface string tension for this value of /3 is cr = 0.0026083 and the dimensional 
reduction scale is exactly 16 lattice spacings [23]. The data correspond to values of the L3 
size ranging from a minimum value of 3/2 the dimensional reduction scale up to 3 times. 

In the table we report, from left to right, the lattice sizes Li.2,3 , the numerical estimate 
-^num for the free energy (with its statistical error), the theoretical estimate F according to 
eq. (4.1), the estimate -f^(o,o) obtained in the full dimensional reduction limit by truncating 
eq. (4.1) to A; = A;' = 0, which corresponds to our 2d formula eq. (2.17), and the estimate 
Fist in which also the first excited states, + A;' = 1, are kept into account. Finally, as 
a comparison, in the last two columns we report the one loop and two loop perturbative 
approximations to the whole NG action discussed, for instance, in [9, 10]. 

In Tab. 2 we report the same set of data evaluated at /3 = 0.226102, corresponding to 
a = 0.0105254 and to a dimensional reduction scale of 8 lattice spacings [23]. 

Looking at these tables one can see that the 2d expression -^(0,0) always gives a very 
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good approximation of the whole NG result F: the difference between the two is always 
smaller than the statistical error on the numerical data. One can also notice that the 2d ap- 
proximation is more reliable when the dimensional reduction scale ^/tt/JSo) is smaller, i.e. 
in Tab. 2, and that within each data set it becomes slightly worse as L3 increases. Looking 
at the seventh column we see that the deviation from the full NG result is completely due 
to the first excited state. 

In fact it is easy to evaluate the gap A between the lowest state and the first excitation 
which appears at the exponent of the asymptotic expansion of the Bessel function in eq. 
(4.1). In the large c7L| limit, the first excited states correspond to A; = 1, /c' = or A; = 0, 
k' = 1 and A takes the value 2itL2/L^. When 7r/3 ^ aL\ ^ tt the first excited state is 
instead the one with A; = 1 , A;' = 1 and the variation of the gap is 

We see that the dominating effect is due to the asymmetry ratio L2 / -L3 (which in our data 
ranges from 8/3 to 4/3) and that a subdominant role is played by the aL^ combination, in 
complete agreement with the results reported in Table 1 and 2. 

Given this agreement between 2d and whole NG results it is very interesting to compare 
them both to the Montecarlo data. The agreement is rather good for large values of F, see 
in particular the results reported in Tab. 2. As F decreases, however, the NG predictions 
compare less favorably to the data, see Table 1. As a matter of fact, in this regime the 
Monte Carlo data seem to better agree with the two loop perturbative expansion than with 
the whole NG result (and its 2d approximation). This fact has already been discussed in 
[16] and is most probably due to the fact that the two loop result is actually consistent also 
at the quantum level, in the sense that at this order the contribution due to the Liouville 
field vanishes [24, 25, 26, 27]. 

As expected, these problems disappear if one works exactly in two dimensions, and in 
fact the direct calculation of the interface partition function in the 2d Ising model reported 
in sec. 3 shows that the expression eq. (2.17) obtained in section 2 should be valid to all 
orders. 



6. Conclusions 

In this paper, we addressed the description of some universal properties of interfaces in 
two-dimensional physical systems. Using a simple model in which one assumes an action 
proportional to the length of the interface, in analogy with the 3d capillary wave model, 
we proposed a general expression for the interface free energy. This expression agrees with 
an exact calculation in the case of the 2d Ising model. We also shew that our 2d result 
represents the dimensional reduction of the interface 3d free energy obtained using the 
capillary wave model, which is tantamount to the Nambu-Goto effective string description, 
when the effects of the conformal anomaly are neglected. Finally, we compared the 3d 
and 2d theoretical predictions to Monte Carlo numerical estimates for the interface free 
energy in the case of the 3d Ising model. In the range of values that we studied the 
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difference between the 3d and 2d estimates is always smaller than the statistical error of 
the numerical data. Moreover we were able to show that this difference is completely due 
to the first excited state of the Nambu-Goto action. 
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